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Based on the phase-space information obtained from a multi-phase transport model within the 
string melting scenario for strange and antistrange quarks, we study the yields and transverse 
momentum spectra of (j) mesons and Q {Q~ + f2+) baryons as well as their anisotropic flows in 
Au-|-Au collisions at RHfC using a dynamical quark coalescence model that includes the effect due 
to quark phase-space distributions inside hadrons. With current quark masses and fixing the (j) and 
n radii from fitting measured yields, we first study the ratio of the yield of baryons to that of <j!> 
mesons as well as their elliptic and fourth-order fiows as functions of their transverse momentum. 
How the elliptic and fourth-order fiows of (j> mesons and fl baryons are related to those of strange 
and antistrange quarks is then examined. The dependence of above results on (f> and 51 radii as well 
as on the strange quark mass is also studied. 

PACS numbers: 12.38.Mh, 5.75.Ld, 25.75.-q, 24.10. Lx 



I. INTRODUCTION 

Recently, there is a lot of interest in using the quark 
coalescence or recombination model to understand the 
experimental data from heavy-ion collisions at the Rel- 
ativistic Heavy Ion Collider (RHIC). As shown in Refs. 
P, 0, ^ iJi] , the quark coalescence model can explain suc- 
cessfully the observed anomalously large enhancement of 
baryon to meson ratio at intermediate transverse mo- 
menta and scaling of the elliptic flow of identified hadrons 
according to their valence quark numbers. Most of these 
applications of the quark coalescence model were based 
on a simple momentum-space coalescence in which only 
quarks with same momentum can coalesce into hadrons. 
Also, the phase-space information of quarks in the par- 
tonic matter produced in relativistic heavy-ion collisions 
is usually taken from a schematic fireball model. Al- 
though the fireball model was still used in Ref. (l| , the ef- 
fect due to quark momentum spread and spatial distribu- 
tion inside hadrons was included in calculating the prob- 
abilities for quarks to coalesce to hadrons. Only in Ref. 
1^ was the coalescence probability calculated also with a 
quark phase-space information from a dynamical parton 
cascade model. Results from such a dynamical quark 
coalescence model indicate that the phase-space struc- 
ture of quarks at freeze out plays an important role in 
hadron production from the quark-gluon plasma. These 
studies were mainly concerned with hadrons which con- 
sist of light quarks and/or heavy charm quarks. In the 
present study, we use the dynamical quark coalescence 
model, that is based on the quark phase-space informa- 
tion from a multi-phase transport (AMPT) model within 
the string melting scenario and includes the quark struc- 
ture of hadrons, to study instead the production and 
anisotropic flow of (j> mesons and (O^ + Q,^) baryons 
that consist of strange quarks in Au+Au collisions at 
RHIC. 



The study of (j) meson and Q baryon production in rel- 
ativistic heavy-ion collisions is a topic of great interest as 
one of the signatures for the quark-gluon plasma (QGP) 
produced in relativistic heavy-ion collisions is enhanced 
production of hadrons consisting of strange and/or antis- 
trange quarks. The enhancement occurs because masses 
of strange quarks are comparable to the temperature of 
the QGP and are thus expected to be abundantly pro- 
duced from quark and gluon inelastic scattering once the 
QGP is formed in the collisions 0, Q- Since the </> me- 
son carries hidden strangeness (ss) and the f2 baryon 
consists of three valence strange quarks (sss) or antis- 
trange quarks (sss), their production in heavy-ion colli- 
sions would also be enhanced if the QGP is formed. 

Furthermore, anisotropic flows of 4> mesons and 
baryons in relativistic heavy-ion collisions are useful for 
understanding the collective dynamics of strange quarks 
in the produced QGP H, Q ■ R is known that anisotropic 
flows are sensitive to the properties of dense matter 
formed du ring the early stage of heavy-ion collisions 
Id .liLIT^IlSlI This sensitivity not only exists in the el- 
hptic flow uTEllIllilllllllllillSlIS but also in the 
smaller higher-order anisotropic flows, such as the fourth- 
order anisotropic flow V4 22, 23, 24, 25, 26, 27, 28, 29^. 
Moreover, observed hadron elliptic flows in heavy-ion col- 
lisions at RHIC were found to satisfy the valence quark 
number scaling, i.e., the elliptic flow per quark is the same 
at same transverse momentum per quark. As shown in 
Refs. HSi, such a scaling of hadron elliptic flows 
according to their valence quark numbers can be un- 
derstood in the quark recombination/coalescence model. 
Other scaling relations among hadron anisotropic flows, 
such as V4{pt) ^ VniPT) , have also been observed in ex- 
perimental data |26l l27l | , and they have been shown in 
the quark coalescence model to relate to similar scaling 
relations among quark anisotropic flows [2^ [l^. Be- 
cause of their small scattering cross sections with other 
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hadrons HHHIlilsllllllllll, mesons and fl 
baryons are not much affected by rescattering effects in 
later hadronic stage of the colhsion and are thus expected 
to provide more direct information on the properties of 
the QGP and how they are produced during hadroniza- 
tion iHllSlllili^. 

Our study shows that in heavy-ion coUisions at RHIC 
the ratio of the yield of baryons to that of cj) mesons at 
intermediate transverse momenta is strongly enhanced 
relative to that at low transverse momenta, similar to 
that observed in the data for the proton to pion ratio. 
Also, the elliptic flows V2{pt) of 4> mesons and fl baryons 
follow approximately the valence quark number scaling. 
Their valence quark number scaled elliptic flows deviate, 
however, strongly from the elliptic flow of strange and 
antistrange quarks. For the forth-order anisotropic flow 
V4{pt), it scales with the square of the elliptic flow for 
both (j) mesons and baryons. We have also studied 
the dependence of above results on the radii of 4> meson 
and baryon as well as on the strange quark mass. For 
physically reasonable radii, both yields of 4> mesons and 
f2 baryons increase with increasing radii. Their yields 
eventually decrease as the radii become unrealistically 
large. The scaled elliptic flows of both (j> mesons and il 
baryons increase with increasing radii and approach that 
of strange and antistrange quarks for very large radii as 
in the naive momentum-space quark coalescence model. 
Changing strange and antistrange quark masses to their 
constituent quark masses affects significantly the yields of 
(j) mesons and 51 baryons but not much their anisotropic 
flows. Our results therefore suggest that both the dynam- 
ical phase-space information of partons at freeze out and 
quark structure of hadrons are important for (f) meson and 
57 baryon production in relativistic heavy-ion collisions. 
We note that our study based on the dynamical quark 
coalescence model for (j) meson production and elliptic 
flow in relativistic heavy-ion collisions is different from 
that of Ref . Q , where </> mesons were reconstructed from 
the invariant mass distributions of kaons and antikaons 
from the AMPT model. 

The paper is organized as follows. In Section ^ we 
briefly review the AMPT model and present detailed in- 
formation on the dynamics of quarks and antiquarks (also 
called partons) in heavy-ion collisions at RHIC. In Sec- 
tion mil the dynamical quark coalescence model based 
on parton phase-space distributions at freeze out is de- 
scribed for (j) meson and Q baryon production. This also 
includes the construction of the quark Wigner phase- 
space functions inside the (j) meson and 51 baryon. The 
dynamical quark coalescence model is then used in Sec- 
tion IIVI to study the yields and transverse momentum 
spectra as well as the anisotropic flows of (j) mesons and 
n baryons in Au-|-Au collisions at RHIC. In Section Ivl 
dependence of above results on the radii of (j) meson and 
57 baryon is studied. The effect due to change in the 
strange quark mass is studied in Section IVII Finally, we 
conclude with a summary in Section IVIII 



II. PARTONIC DYNAMICS IN HEAVY-ION 
COLLISIONS AT RHIC 

The dynamical quark coalescence model used in 
present study requires the space-time and momentum 
information of quarks and antiquarks during hadroniza- 
tion of the produced QGP in relativistic heavy-ion col- 
lisions. In particular, the dynamics of strange and anti- 
strange quarks are needed for describing the production 
of (j) mesons and 57 baryons as well as their anisotropic 
fiows, and they are taken from the AMPT model within 
the string melting scenario 



A. The AMPT model 

The AMPT model is a hybrid model that uses mini- 
jet partons from hard processes and strings from soft 
processes in the Heavy Ion Jet Interaction Generator 
(HIJING) model as the initial conditions for mod- 
eling heavy-ion collisions at ultra-relativistic energies. 
Time evolution of resulting minijet partons, which are 
mostly gluons, is described by Zhang's parton cascade 
(ZPC) '43| model. At present, this model includes only 
parton-parton elastic scatterings with in-medium cross 
sections derived from the lowest-order Born diagrams 
and having magnitude and angular distribution fixed 
by treating the gluon screening mass as a parameter. 
After minijet partons stop interacting, they are com- 
bined with their parent strings, as in the HIJING model 
with jet quenching, to fragment into hadrons using the 
Lund string fragmentation model as implemented in the 
PYTHIA program . The final-state hadronic scatter- 
in^ s are modeled by a relativistic transport (ART) model 

Since the initial energy density in heavy-ion collisions 
at RHIC is much larger than the critical energy density at 
which the hadronic matter to quark-gluon plasma transi- 
tion would occur |^ ,47^ , the AMPT model has been 
extended by converting initial excited strings into par- 
tons In this string melting scenario, hadrons that 
would have been produced from string fragmentation are 
converted instead to their valence quarks and/or anti- 
quarks. Interactions among these quarks are again de- 
scribed by the ZPC parton cascade model. Since inelastic 
scatterings are at present not included, the resulting par- 
tonic matter consists of only quarks and antiquarks from 
melted strings. To take into account the effect of stronger 
scattering among gluons if they were present, the scat- 
tering cross sections between quarks and antiquarks are 
taken to be the same as those for gluons. These quarks 
and antiquarks are converted to hadrons when they stop 
scattering with other partons. For a parton cross section 
of 10 mb, which is used in present study, this criteria for 
hadronization is not too different from the one based on 
the condition that the local energy density of hadronizing 
partons is at the critical energy density given by lattice 
QCD calculations. The transition from the partonic mat- 
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ter to the hadronic matter is achieved using a simple co- 
alescence model, which combines two nearest quark and 
antiquark into mesons and three nearest quarks or anti- 
quarks into baryons or anti-baryons that are close to the 
invariant mass of these partons. Details of the AMPT 
model can be found in Rcf. j4l| . 

In the present study, we use the AMPT model in the 
string melting scenario with the default current quark 
mass of 9.9 MeV for d quark, 5.6 MeV for u quark, and 
199 MeV for s quark and a constant parton scattering 
cross section of 10 mb. We note that using parton scat- 
tering cross sections of 6-10 mb, the AMPT model with 
string melting was able to reproduce both the centrality 
and transverse momentum (below 2 GeV/c) dependence 
of the elliptic flow 18] and pion interferomctry 49] mea- 
sured in Au-f Au collisions at ^/snn = 130 GeV at RHIC 
[53,1^ as well as measured pr dependence of both V2 and 
V4 of mid-rapidity charged hadrons in the same collision 
at y/sjn^ = 200 GeV p^ . 

B. Spectrum of strange and antistrange quarks 
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FIG. 2: (Color online) Transverse momentum distributions of 
mid-rapidity strange and antistrange quarks (s-|-s) as well as 
light quarks and antiquarks {u + u + d + d) at freeze out in 
Au+Au collisions at yijvjv ~ 200 GeV and 6 = 8 fm.. 
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FIG. 1: (Color online) Rapidity distributions of strange and 
antistrange quarks (s -I- s) as well as light quarks and anti- 
quarks {u + u + d + d) at freeze out in Au+Au collisions at 
V'si^ = 200 GeV and fe = 8 fm. 

We first show in Fig. ^ the rapidity distributions of 
strange quarks and antiquarks (s + s) as well as light 
quarks and antiquarks {u + u + d + d) at freeze out, i.e., 
stop scattering with other partons, in Au+Au collisions 
at y^SNN = 200 GeV and 6 = 8 fm. It is seen that the 
rapidity density dN/dy at mid-rapidity is about 94 for 
strange quarks and antiquarks and about 770 for light 
quarks and antiquarks. Since inelastic partonic scatter- 
ings are absent in the present model, all these quarks and 
antiquarks are from melted strings during initial stage of 
the collisions. 



The transverse momentum distributions of midrapid- 
ity strange and antistrange quarks as well as light quarks 
and antiquarks at freeze out in the same collision are 
shown in Fig. |21 It clearly shows a quark mass effect as 
strange and antistrange quarks display a stiffer transverse 
momentum distribution at low pr compared to that of 
light quarks and antiquarks. Fitting the transverse mo- 
mentum spectrum by a Boltzmann distribution, we find 
that the effective temperature of the partonic matter at 
freeze out is about 161 MeV, which is consistent with the 
predicted critical temperature Tc ~ 150-180 MeV from 
the lattice QCD for the QGP to hadronic matter phase 
transition 153. ISal. 



C. Space-time structure of strange and antistrange 
quarks at freeze out 

The spatial distribution of mid-rapidity strange and 
antistrange quarks at freeze out in the above collision is 
shown in Fig. |3| with the upper panel displaying their 
number density distribution in the transverse plane and 
the lower panel giving corresponding contour plot. Here, 
the transverse plain refers to the x-y plane with the a;-axis 
pointing to the direction of the impact parameter and the 
j/-axis perpendicular to the z-axis as well as the beam 
direction. It is seen that the number density of strange 
and antistrange quarks at freeze out peaks approximately 
at a circle with radius of about 4 fm in the transverse 
plane, implying that these partons mainly freeze out from 
the surface of an expanding fireball where the number 
density is roughly between 1.2 fm~^ and 1.4 fm^^ in the 
transverse plane. 

The space-time structure of z-t correlations for mid- 
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FIG. 3; (Color online) Number density distribution in the 
transverse plane for mid-rapidity strange and antistrange 
quarks at freeze out in Au+Au collisions at y^ijvjv = 200 
GeV and 6 = 8 fm (upper panel) and corresponding contour 
plot (lower panel). 
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FIG. 4: (Color online) Space-time structure of mid-rapidity 
strange and antistrange quarks at freeze out in Au-|-Au colli- 
sions at ^/snn = 200 GeV and fe = 8 fm. 
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rapidity strange and antistrange quarks at freeze out in 
this collision is shown in Fig. 21 Instead of a sudden 
freeze out, these partons are seen to freeze out contin- 
uously in the proper time r, with most between t — A 
and 6 fm/c. This can be more clearly seen from their 
freeze out rate as a function of proper time r or the time 
t in the nucleus-nucleus center-of-mass system as shown 
in Fig. O It shows that the freeze out rate of midrapidity 
strange and antistrange quarks peaks at r = 5 fm/c or 
t = 5.5 fm/c. 

The AMPT model also provides information on the 
time evolution of parton energy density. This is shown 
in Fig. for the central cell of the partonic matter in 
above collision. As in Refs. 0,E3i the central cell is 
taken to have a transverse radius of 1 fm and a longitu- 



FIG. 5: (Color online) Freeze out rate of mid-rapidity strange 
and antistrange quarks as a function of proper time r as well 
as the time t in the nucleus-nucleus center-of-mass system. 



dinal dimension of 5% of the time t after the two nuclei 
have fully overlapped in the longitudinal direction. It is 
seen that the central energy density is about 0.6 GeV/fm'^ 
a,t t — 5.5 fm/c when most of partons freeze out. This 
energy density is consistent with the critical energy den- 
sity of about 0.66 GcV/fm"^ predicted by the lattice QCD 
for the quark-gluon plasma to hadronic matter transition 

MM- 
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FIG. 6: (Color online) Time evolution of the energy density 
in the central cell of partonic matter in Au+Au collisions at 
^^sWn = 200 GeV and fe = 8 fm. 

D. Anisotropic flow of strange and antistrange 
quarks 

Anisotropy in the transverse momentum distribution 
of particles in non-central heavy-ion colhsions is gener- 
ated by the pressure anisotropy in the initial compressed 
matter ^4, 55] and thus depends on the geometry and en- 
ergy density as well as the properties of produced matter 
during the early stage of these collisions. For (j) mesons 
and Q baryons, their anisotropic flows are essentially de- 
termined by those of strange and antistrange quarks at 
freeze out and their hadronization mechanism. As dis- 
cussed previously, the effect due to later hadronic scat- 
tering is small because the interactions of mesons and 
O baryons in the hadronic matter are relatively weak 
and also the pressure becomes more isotropic during the 
hadronic stage of the collisions. 

From the momentum distribution of strange and antis- 
trange quarks at freeze out, their elliptic and fourth-order 
anisotropic flows can be evaluated according to following 
single-particle averages: 

Mpt) = M~'P'fi+A (2) 
\ Pt / 

where Px and Py are, respectively, projections of their 
momentum in and perpendicular to the reaction plane, 
defined by the beam and impact parameter axes, and 
Pt = [Px + P^y^^ is the transverse momentum. 

In Fig. 13 we show the pr dependence of V2 and 
of mid-rapidity strange and antistrange quarks (s -I- s) at 
freeze out in Au-|-Au collisions at ^/snn = 200 GeV and 
b = 8 fm. It is seen that strange and antistrange quarks 
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FIG. 7: (Color online) Transverse momentum dependence 
of elliptic flow U2 (solid squares) and forth-order anisotropic 
flow D4 (open squares) of mid-rapidity strange and antis- 
trange quarks (s + s) at freeze out in Au+Au collisions at 
y/sNN ~ 200 GeV and b — 8 fm. The thick solid line rep- 
resents 0.85«2 of strange and antistrange quarks, and open 
triangles are V2 of mid-rapidity light up and anti-up quarks 
(it -f It). 

exhibit not only a strong V2 but also a non-negligible V4 
that approximately scales with with a scaling coef- 
ficient of 0.85, i.e, V4 ~ 0.85w| as shown by the thick 
solid line in Fig. [7| Compared to the V2 of mid-rapidity 
light up and anti-up quarks (u + u), shown by triangles 
in Fig. \7\ V2 of heavier strange and antistrange quarks 
has a smaller value at low px but a larger value at high 
Pt, similar to the mass ordering of hadron elliptic fiows 
in the hydrodynamic model which assumes that the pro- 
duced matter is in local thermal equilibrium and thus 
develops a large collective radial flow. However, instead 
of a continuing increase of V2 with respect to pr as in 
the hydrodynamic model, V2 of partons in the transport 
model saturates at a maximum value when their pT be- 
comes large, indicating that high momentum partons do 
not reach thermal equilibrium with the bulk of the par- 
tonic matter. 



III. A DYNAMICAL QUARK COALESCENCE 
MODEL FOR AND PRODUCTION 

In the coalescence model, the probability for forming 
a bound cluster from a many-particle system is deter- 
mined by the overlap of the wave functions of coalescing 
particles with the internal wave function of the cluster. 
Its validity is based on the assumption that coalescing 
particles are statistically independent and the binding 
energy of formed cluster as well as the quantum dynam- 
ics of the coalescing process play only minor roles [56l |. 
In the present study, we assume that correlations among 
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partons at freeze out are weak and binding energies of 
formed hadrons can be neglected. Furthermore, the coa- 
lescence model is considered as a perturbative approach, 
valid only if the number of partons coalesced into hadrons 
is small compared with the total number of partons in the 
system. This condition is satisfied for cj) mesons and 51 
baryons produced in relativistic heavy-ion collisions as 
their numbers measured in experiments are indeed sig- 
nificantly smaller than that of kaons and thus that of 
strange and antistrange quarks. 



A. The dynamical coalescence model 

In the dynamical quark coalescence model, the proba- 
bility for producing a hadron from partons in the QGP is 
given by the overlap of parton phase-space distributions 
at freeze out with the parton Wigner phase-space func- 
tion inside the hadron. For a QGP containing A partons, 
the momentum distribution of a M -parton hadron can be 
expressed as [H^ |E3 



dN 



M 



= G 



A 
M 



1 f 



xp 
xS{K 



(ki 



kA/))dridki ■ 



■dvMclkM, (3) 



where fi is the parton phase-space distribution func- 
tions at freeze-out; is the M-parton Wigner phase- 
space function inside the hadron; ri^,--- ,rijj_j and 
kij, • • • ,kijjj_j are, respectively, the M — 1 relative co- 
ordinates and momenta taken at equal time in the M- 
parton rest frame (K = 0); and G is the statistical factor 
for the M partons to form the hadron. 

In transport model simulations of heavy- ion collisions, 
the multiplicity of a M-parton hadron produced from 
the dynamical quark coalescence model is then given by 

El El 



N 



M 



G 



E 



■ dr d\c 



■il>i2>--->4Af 



-1 J ^iM-l )) 7 



(4) 



where (• • • ) denotes event averaging and the sum runs 
over all possible combinations of M partons. 

To determine the statistical factor G for the (j) meson, 
we note that its quark wave function in the color-spin- 
isospin space can be expressed as a linear combination 
of all possible orthogonal flavor, color, and spin basis 
states [l|, . The probability for a strange quark s and 
an anti-strange quark s to form a hadron with quan- 
tum numbers corresponding to a (f) meson with the z- 
component of its spin equal to 1 is simply given by the 
probability of finding the two s and s in any one of these 
color-spin-isospin basis states, i.e., 1/3^ x 1/2^ = 1/36. 
Including also the possibility of forming a (j) meson with 
the z-component of its spin equal to —1 and triples 



the probability. As a result, the statistical factor G for 
the (j> meson is 1/12. A similar consideration leads to a 
statistical factor G = 1/54 for the H.^ and f2+ baryons. 



B. Parton phase-space distributions 

The M-parton phase-space distribution function 
nfii/(ri,ki) in Eq. Q refers to partons with spatial 
coordinates and momenta at equal time in the rest frame 
of the M-parton cluster. Since parton momenta in the 
AMPT model are given in nucleus-nucleus center-of-mass 
system, a Lorentz transformation is performed to obtain 
their momenta in the rest frame of the M-parton clus- 
ter. Furthermore, partons freeze out at different times as 
shown in Figs. 01and[5| The same Lorentz transformation 
is thus used in obtaining their space-time coordinates in 
the rest frame of the M-parton cluster. To determine the 
spatial coordinates of these partons at equal time in their 
rest frame, partons in the cluster that freeze out earlier 
are allowed to propagate freely, i.e., with constant veloc- 
ities given by the ratio of their momenta and energies in 
the rest frame of the cluster, until the time when the last 
parton in the cluster freezes out. 



C. Quark Wigner phase-space functions inside (f> 
and 

To determine the quark Wigner phase-space functions 
inside hadrons requires knowledge of their quark wave 
functions. For the (j) meson, we take its quark wave func- 
tion to be that of a spherical harmonic oscillator, i.e.. 



V'(ri,r2) = l/(7ra^)3/4exp(-rV(2a^)), 



(5) 



in terms of the relative coordinate r = ri — r2 and the size 
parameter . The above normalized wave function leads 

to a root mean-square radius — (r^)^^^ = (3/8)^/^(70 
for the meson. 

The quark Wigner phase-space function inside the 
meson is obtained from its quark wave function by 



r,k) 



1/; -I- 5: j ^* f^r - y ) exp(-ik • R)d3R 



2;„2 



Sexp 2 - ajz 



(6) 



where k = (ki — k2)/2 is the relative momentum between 
s and s quarks. 

For rj" and fl^ baryons, their quark wave functions are 
taken to be the same and are given by that of a spherical 
harmonic oscillator as well |EalH3i i-^-; 



V;(ri,r2,r3) = (3^2^^)-3/4exp 



(7) 
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in terms of the relative coordinates p and A , and the 
size parameter ctq. Here we have used the usual Jacobi 
coordinates for a three-particle system i.e., 




Ill 

%/6 V6 %/6 




(8) 



where R is the center-of-mass coordinate of the three 
quarks or antiquarks. 

Using dridv2dY^ — S'^^^dUdpdX, it is easy to check that 
the wave function given by Eq. ^ is normalized to one. 
From the relation (ri — R) + (r2 — R) + (fa — R) = 

+ A'^, the root mean-square radius Rn of the fl baryon 
is given by 



Rn = 



|V'(ri,r2,r3)p33/2dpdA 



-,1/2 



(9) 

The quark Wigner phase-space function inside the 
baryon is obtained from its quark wave function via 



Pn (P, A,kp,kA) 

/ 1 
tp p 



~ » cxp 



— ,A+— (^P- — ,A- — 

exp(-ikp • Ri) exp(-ikA • R2)3^/^dRi<iR2 
P^ + A^- 



cxp(-(k^ + kl)a^), (10) 



where kp and k^ are relative momenta, which together 
with the total momentum K are defined by [s^ 



of 0-5%. Using the above described dynamical quark co- 
alescence model, we find that the STAR (PHENIX) data 
can be approximately reproduced with a meson size 
parameter ~ 1.06 (0.77) fm, which gives a reasonable 
(j) meson root-mean-square radius of about — 0.65 
(0.47) fm. For the impact parameter 6 = 8 fm as mainly 
considered in present study, the AMPT model gives the 
number of participant nucleons Npart ~ 160 which cor- 
responds to a centrality of about 31%. The resulting 
meson rapidity density dN/dy at mid-rapidity obtained 
with the (j) meson radius R^ = 0.65 (0.47) fm is about 3.3 
(2.0), which is also in good agreement with experimental 
data (3&i39l|. 

For Q baryons, there are only preliminary experimen- 
tal data on their dN/dy at mid-rapidity in central colli- 
sions of Au-I-Au at ^/sNN — 200 GeV, and its value is 
about 0.64 Fitting the result from the dynamical 

quark coalescence model for the same collision at 6 = 
fm leads to an baryon size parameter crn « 1.2 fm, 
corresponding to an Q baryon root-mean-square radius 
of about Rq « 1.2 fm, which is somewhat large but not 
unrealistic. For more peripheral collisions at impact pa- 
rameter 5 = 8 fm, an SI radius of Rq — 1.2 fm leads to a 
rapidity density dN/dy of about 0.26 at mid-rapidity. 

In the following, we use the above determined (j) and f2 
root mean-square radii, i.e., R^ = 0.65 fm or 0.47 fm for 
the (j) meson while Rq = 1.2 fm for the O baryon unless 
stated otherwise. 



IV. (I) AND n PRODUCTION IN AU+AU 
COLLISIONS AT = 200 GeV 




K 



with ki, k2, and ka being the momenta of the three 
quarks. 
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Using the phase-space information and the dynamical 
coalescence model described above, we can now study the 
transverse momentum spectra of (j) mesons and fl baryons 
as well as their anisotropic fiows in relativistic heavy-ion 
collisions. 



D. (j) and size parameters 

The two parameters and an in the quark Wigner 
phase-space functions inside the meson and Q baryon 
are related to their root-mean-square radii. Since the lat- 
ter are not known empirically, we take them as adjustable 
parameters and fix them by fitting measured yields of (j) 
mesons and SI baryons in relativistic heavy- ion collisions. 

For (f) mesons, experimental data on their rapidity den- 
sity dN/ dy at mid-rapidity and its centrality dependence 
are available for Au+Au collisions at ^/snn — 200 GeV 
[38L Isgj. There is, however, a significant difference be- 
tween the STAR and PHENIX data for the meson yield. 
For instance, the </> meson dN/dy at mid-rapidity is about 
7.7 in the STAR data while 4.5 in the PHENIX data for 
Au-|-Au collisions at ^sjvAf = 200 GeV with a centrality 



A. Transverse momentum spectra of and Q, 

Fig. |S1 shows the transverse momentum spectra of 
midrapidity 4> mesons and SI baryons as well as their ratio 
in Au-|-Au collisions at ^s^n — 200 GeV with 6 = fm 
(left panels) and 6 = 8 fm (right panels), and R^ — 0.65 
fm or 0.47 fm and Rq = 1.2 fm. For comparison, we also 
include in upper left panel of Fig. |H1 the experimental 
data for (j> mesons in 0-5% centrality bin from the STAR 
collaboration [s^l (solid stars) and in 0-10% centrality bin 
from the PHENIX collaboration 39] (open stars) as well 
as the preliminary data for S7 baryons in 0-10% central- 
ity bin from the STAR collaboration |0 (solid triangles) . 
Also shown in lower left panel by open triangles are data 
on the fl/cj) ratio in 0-10% centrality bin, obtained from 
a linear interpolation of the (p meson data in upper left 
panel. 
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FIG. 8; (Color online) Transverse momentum dependence 
of midrapidity (j> mesons and fl (f2~ + Q'^) baryons (upper 
panels) and their ratio (lower panels) in Au-f An collisions at 
-ysivjv ~ 200 GeV with 6 = fm (left panels) and 6 = 8 fm 
(right panels). Data for (j> mesons in -5% centrality bin from 
the STAR collaboration [s^ l (solid stars) and in 0-10% cen- 
trality bin from the PHENIX collaboration [soll (open stars) 
as well as preliminary data for fl baryons in 0-10% central- 
ity bin from the STAR collaboration (solid triangles) are 
shown in upper left panel, while those for the fl/cp ratio in 
0-10% centrality bin are shown in lower-left panel. 




flows are given by ratios of hadron transverse momentum 
spectra, predictions of anisotropic flows as well as the ra- 
tio of hadron transverse momentum spectra frorn present 
AMPT model are expected to be more reliable [Sill. 

The ratio of midrapidity Q baryons to (j) mesons shown 
in the lower panels of Fig. |Slis seen to increase apprecia- 
bly from low pt to intermediate pr of about 2.5 GeV/c. 
The enhancement factor is about 7 and 4 for = fm 
and b — 8 fm, respectively. This enhancement is a result 
of parton coalescence as in the case of observed anoma- 
lously large anti-proton to pion ratio of about 0.8 and 0.4 
^■t Pt ^ 2.5 GcV/c in central and mid-peripheral (cen- 
trality of about 30%) Au-|-Au collisions at ^/sNN = 200 
GeV, respectively [63. As explained in Refs. 0,0>I11j 
in the coalescence model baryons with transverse mo- 
mentum Pt are mainly formed from quarks with trans- 
verse momenta ^ Pt/'^i while mesons with same trans- 
verse momentum are mainly produced from partons with 
transverse momenta ^ pt/'^- Since quark transverse mo- 
mentum spectra decrease with pT, it is more favorable to 
produce high transverse momentum baryons than mesons 
if hadrons are produced from the QGP through the co- 
alescence of quarks. However, the r2+/(/) ratio is only 
about 0.23 and 0.17 at pt « 2.5 GeV/c for & = fm and 
b — 8 fm, respectively, and is significantly smaller than 
the anti-proton to pion ratio observed in experiments. 
For b — fm, our result with — 0.65 fm reproduces 
reasonably the measured fl/cj) ratio as shown in lower left 
panel of Fig. |H1 



It is seen from upper left panel of Fig. |S1 that the 
calculated transverse momentum spectrum of baryons 
is stiffer than that of (j) mesons as expected from their 
different masses. Both calculated transverse momentum 
spectra of mesons and baryons are, however, sig- 
nificantly softer than the data. With = 0.65 (0.47) 
fm, the mean transverse momentum calculated from the 
(j) meson transverse momentum spectrum is about 0.72 
(0.73) and 0.68 (0.70) GeV/c for & = fm and = 8 fm, 
respectively, compared with the experimental value of 
about 0.85-1.1 GeV/c ||3^|3^. For SI baryons, their mean 
transverse momentum is about 0.93 and 0.83 GeV/c for 
& = fm and b = 8 fm, respectively, and are larger than 
those of (j) mesons. 

We note that hadron transverse momentum spectra 
obtained from the AMPT model with string melting are 
generally softer than those measured in experiments j61j| . 
Since hadron spectra refiect those of partons in the par- 
ton coalescence model, soft hadron spectra are due to soft 
parto n sp ectra during the partonic stage. As discussed in 
Ref. (63 , the latter is a result of the small current quark 
masses used in the parton cascade of the AMPT model 
with string melting, which make their transverse momen- 
tum spectra less affected by radial collective flow than 
hadrons in the default AMPT model. Since anisotropic 
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FIG. 9: (Color online) Transverse momentum dependence of 
the ratio of midrapidity (f} meson and baryon yields in cen- 
tral (& = fm) to that in mid-peripheral {b — 8 fm) Au-j-Au 
collisions at ^snn — 200 GeV. 

In Fig. we show the pt dependence of the ratio of 
midrapidity (j) meson or f2 baryon yield in central (& = 
fm) to corresponding one in mid-peripheral (6 = 8 fm) 
collisions of Au-f Au at ^/snn 



200 GeV. It is seen 
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that the pt dependence is very different for (f> mesons 
and Q baryons with the latter having a smaller value at 
low Pt while a larger value at high px than (j) mesons. 
Our results thus indicate that nuclear suppression at in- 
termediate momenta is weaker for baryons than (p 
mesons. This feature is consistent with the measured 
Pt dependence of so-called Rcp (the ratio of binary- 
collision-scaled Pt spectra in central and peripheral col- 
lisions) for p + p and pions [6^ . 



B. Anisotropic flows of <j) and fl 

In heavy- ion collisions at RHIC, experimental data 
have indicated that there exist not only strong ellip- 
tic flow but also a clear fourth-order anisotropic flow 
V4 for charged hadrons [2^ l2^ . An interesting find- 
ing about W4 is that it scales with the square of V2, i.e., 
'Vi{PT) ~ ''^il-PT)- III tli6 present subsection, both ellip- 
tic and fourth-order anisotropic flows of mesons and 
VL baryons are studied to see if they also show similar 
features. 
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FIG. 10: (Color online) Transverse momentum dependence 
of anisotropic flows 112 and va of mid-rapidity <j) mesons (left 
panel) and Q, baryons (right panel) produced in Au-|-Au col- 
lisions at ^snn = 200 GeV and b = 8 fm. Solid and dashed 
lines in left panel are, respectively, l.lvi and 1.2v2 for (f> 
mesons, while solid line in right panel is 0.7v2 for fl baryons. 

In left panel of Fig. ^1 we show the pt dependence 
of anisotropic flows V2 and Vi of mid-rapidity (j) mesons 
produced in Au-|-Au collisions at ^^snn = 200 GeV and 
6 = 8 fm. Solid and open squares are, respectively, the 
meson V2 and V4 for a (f) meson radius = 0.65 fm, while 
solid and open triangles are, respectively, those for = 
0.47 fm. It is seen that the scaling relation W4(pt) ^ 
V2{pt) is satisfied in both cases as shown by the solid 
(l.lw|) and dashed (1.2i;|) lines in the same figure. It is 
interesting to see that the scaling coefficient of 1.1 or 1.2 
is the same as that extracted from the data for charged 



hadrons |26l |27| . Similar results for fl baryons using the 
radius Rq = 1.2 fm are shown in right panel of Fig. [TUl 
The scaling relation v^Ipt) ^ ''^il-P'r) is again satisfied 
for f2 baryons but with a smaller scaling coefficient than 
that for ip mesons as shown by the solid line (0.7w|) in 
the right panel. 

Based on the naive momentum-space quark coales- 
cence model that only allows quarks with equal momen- 
tum to form hadrons 0, this scaling relation between 
hadron V2(j)t) and V4{pt) has been attributed to a similar 
scaling relation between those of quarks 28, 29] . Neglect- 
ing the small contribution of higher-order anisotropic 
flows (higher than forth-order), this model gives for mid- 
rapidity hadrons 



^^2,Af(2PT) 

V4,Bi3pT) 
^iBi^PT) 



1 1 V4,q{PT 

4 



2i;I4pt)' 




(12) 



where Vn^M{pT), Vn^sipr) and Vn,q{pT) denote, respec- 
tively, the meson, baryon, and quark anisotropic flows. 
The hadron anisotropic flows thus satisfy the scaling re- 
lation Vi{pT) ^ 'V2{pt) if a similar scaling relation is 
satisfied by quark anisotropic flows. 

We have already shown in Fig. [3 that the scaling rela- 
tion ViipT) ~ "f^KPr) is satisfied by mid-rapidity strange 
and antistrange quarks (s-l-s) with a scaling coefficient of 
about 0.85, i.e., U4 

,qiPT)/vlqiPT) ~ 0.85. Eq. (O then 
leads to a scaling coefficient of V4^^^{pt)Iv2 ^{pt) ~ 0.68 
for (j) mesons and Vi^a{PT)/v2 q{pt) ~ 0.62 for fJ baryons. 
Comparing with the scaling coefficients of 1.1 or 1.2 for 
(j) mesons and 0.7 for f2 baryons from the dynamical 
quark coalescence model, the predicted value from the 
naive momentum-space quark coalescence model for the 
(j) meson scaling coefficient is significantly smaller. How- 
ever, the dynamical quark coalescence model does give a 
larger scaling coefficient for cf) mesons than for f2 baryons 
as expected from the naive momentum-space coalescence 
model. 

Another interesting and important finding in heavy-ion 
collisions at RHIC is the valence quark number scaling of 
the elliptic flow of identified hadrons, i.e., the elliptic flow 
per valence quark in a hadron is same at same transverse 
momentum per valence quark, i.e., V2,H{PT/nq)/nq = 
V2,q{PT) with V2^H and Uq denoting, respectively, the 
hadron V2 and the number of valence quarks or anti- 
quarks in a hadron. In the naive momentum-space quark 
coalescence model, this scaling has been shown to be sat- 
isfied if high-order anisotropic flows are small j§, ; 2^ . 
The scaling is, however, broken in more general quark 
coalescence models that take into account the quark mo- 
mentum distribution , ,63- .64. .65.] and higher parton Fock 
states Ipg in hadrons as well as the effect of resonance 
decays [63 . It is thus of interest to see if the elliptic flows 
of (j) mesons and VL baryons from present dynamical quark 
coalescence model also satisfy the valence quark number 
scaling. 
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FIG. 11: (Color online) Valence quark number scaled ellip- 
tic flow V2/nq as a function of scaled transverse momentum 
Pt /nq for mid-rapidity (j) mesons and fl baryons produced in 
Au-|-Au collisions at ^/sNN = 200 GeV and fe = 8 fm with 
= 0.65 fm (solid squares) and 0.47 fm (open squares) 
as well as Rn = 1.2 fm (solid line). Results for mid-rapidity 
strange and antistrange quarks (s-|-s) at freeze-out are shown 
by solid stars. 



In Fig. 111! we show the valence quark number scaled 
elliptic flow V2 /nq as a function of scaled transverse mo- 
mentum Pt I Uq for mid-rapidity mesons and baryons 
produced in Au-|-Au collisions at ^Jsnn = 200 GeV and 
6 = 8 fm with = 0.47 or 0.65 fm and Rn = 1.2 fm. 
For comparison, we also include the elliptic flow of mid- 
rapidity strange and antistrange quarks [s + s) at freeze- 
out. It is seen that the elliptic flows of both mid-rapidity 
(j) mesons and f2 baryons from the dynamical quark co- 
alescence model satisfy the valence quark number scal- 
ing. The valence quark number scaled elliptic flows of (j) 
mesons and baryons are, however, significantly smaller 
than that of coalescing strange and antistrange quarks. 
This feature is different from the prediction of the naive 
momentum-space quark coalescence model where the va- 
lence quark number scaled V2 of hadrons is equal to the V2 
of coalescing quarks. Our results therefore indicate that 
both the dynamical phase-space information of quarks 
and the quark phase-space distribution inside hadrons 
play important roles in hadron anisotropic flows, and 
cautions are needed in interpreting the elliptic flow of 
partons from that of hadrons using the naive momentum- 
space parton coalescence model. We note that our results 
are consistent with previous studies on V2 of protons and 
pions based on the parton phase-space information from 
Molnar's parton cascade (MFC) model j^. 



V. HADRON SIZE DEPENDENCE OF (j) AND Q. 
YIELDS AND ANISOTROPIC FLOWS 

In above calculations, we have used (p meson and f2 
baryon size parameters that are fixed from fitting mea- 
sured yields in experiments. Ideally, we would like to 
use the empirically and/or theoretically determined root- 
mean-square radii of </> meson and f2 baryon in the dy- 
namical quark coalescence model. In the absence of such 
information, it is of interest to study how above results 
depend on the sizes of meson and f2 baryon. 



A. Hadron size dependence of (j> and Q. yields 
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FIG. 12: (Color online) Hadron size (root-mean-square ra- 
dius) dependence of rapidity density dN/dy of mid-rapidity 
(j) mesons (solid squares) and Q, (f2~ -I- Q,'^) baryons (open 
squares) produced in Au-|-Au collisions at ^/snn ~ 200 GeV 
with b = fm (upper panel) and fe = 8 fm (lower panel). 

We show in Fig. El the rapidity density dN/dy of mid- 
rapidity (j) mesons and f2 baryons produced in Au+Au 
collisions at ^/SNN — 200 GeV with 6 = fm (up- 
per panel) and 6 = 8 fm (lower panel) as functions of 
their root-mean-square radii. It is seen that both (p me- 
son and Q baryon yields are sensitive to their sizes. For 
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physically reasonable radii, both yields increase with in- 
creasing radii but decrease eventually as the radii become 
unrealistically large. For collisions at impact parameter 
6 = fm, the yield exhibits a maximum at i?^ = 2.2 
fm for (f) mesons and at Rn = 2.5 fm for Q, baryons. A 
similar feature is observed in the case of 6 = 8 fm but 
the maximum is at = 1.75 fm for <j) mesons and at 
i?n = 2 f m for 17 baryons. 

B. Hadron size dependence of (j) and Q. elliptic 
flows and the valence quark number scaling 
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FIG. 13: (Color online) Transverse momentum dependence of 
the V2 of mid-rapidity (j) mesons (left panel) and Q. baryons 
(right panel) in Au-(-Au collisions at ^snn = 200 GeV and 
6 = 8 fm with different values of root-mean-square radii R^f, 
and Rn. 



eters cr^ and trn) are very large. As a result, we expect 
that the valence quark number scaled elliptic flows of (p 
mesons and f2 baryons should approach the V2 of coalesc- 
ing strange and antistrange quarks as their sizes increase. 
This is shown in Fig. 1141 where it is seen that the pt /nq- 
dependence of V2lnq for mid-rapidity 4> mesons and f2 
baryons produced in Au-|-Au collisions at ^/snn = 200 
GeV and 6 = 8 fm with R^ — Rn = 4.5 fm is indeed 
similar to the elliptic flow of mid-rapidity strange and 
antistrange quarks (s -I- s) at freeze out. 
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FIG. 14: (Color online) Same as Fig. [TTl but with R4, = 
Rn = 4.5 fm. 



C. Fourth-order flows of (fi and fl 



We first show in left panel of Fig. ^| the pr- 
dependence of the V2 of mid-rapidity mesons in Au+Au 
collisions at ^/snn = 200 GeV and 6 = 8 fm for different 
values of root-mean-square radius i?^. It is clearly seen 
that the V2 of mesons depends strongly on the (f) meson 
size. It becomes larger with increasing R^ but changes 
little when i?^ > 4.5 fm. Similarly, we show in right panel 
of Fig. El the pT-dcpendence of the f 2 of mid-rapidity 

baryons produced in the same collision with different 
values of Rq. Similar to that of cf) mesons, the V2 of 
baryons also increases with increasing size of ft baryon 
and saturates when Rn > 4.5 fm. The saturation of (f> 
meson and baryon V2 at large sizes observed in Fig. 1131 
is essentially due to the fact that the dynamical quark co- 
alescence model approaches the naive momentum-space 
quark coalescence model, i.e, only partons with same mo- 
mentum can coalesce into hadrons, when the hadron size 
is large enough. This can be easily seen from Eqs. © 
and H1U|) as the quark Wigner phase-space function be- 
comes essentially delta functions in the relative momenta 
of coalescing partons when hadron sizes (i.e., the param- 
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(Color online) Same as Fig. llOl but with R^ = Rq = 
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For the scaling relation 174 (pr) ^ ^^Kpt) among 
hadron anisotropic flows, as discussed above and in Refs. 
[2^ [2^, the dynamical quark coalescence model leads 
to a deviation of the scaling coefficient from that of 
the naive momentum-space quark coalescence model as 
a result of finite hadron sizes and nonzero parton rel- 
ative momenta inside hadrons. The deviation is ex- 
pected to disappear for large enough hadron sizes when 
the dynamical quark coalescence model approaches the 
naive momentum-space quark coalescence model. This 
is demonstrated in Fig. 1151 where we show the p^- 
dependence of the anisotropic flows V2 (solid squares) and 
Vi (open squares) of mid-rapidity mesons (left panel) 
and O baryons (right panel) produced in Au+Au colli- 
sions at ^ysNN = 200 GeV and 5 = 8 fm with large 
hadron sizes, i.e., = i?n = 4.5 fm. Also shown by 
solid lines are 0.68i;2 for cj) mesons in left panel and 0.62'i;| 
for baryons in right panel. It is seen that for such large 
hadron sizes the scaling coefficients indeed approach the 
values expected from the naive momentum-space quark 
coalescence model, i.e., about 0.68 for (p mesons and 0.62 
for baryons. We note that a larger hadron size gener- 
ally gives a smaller scaling coefficient in the dynamical 
quark coalescence model. 



VI. QUARK MASS DEPENDENCE OF AND 
TRANSVERSE MOMENTUM SPECTRA AND 
ANISOTROPIC FLOWS 

Above results were obtained with the current strange 
quark mass of 199 MeV, which is the default value in 
the AMPT model with string melting. Usually, the con- 
stituent quark mass (about 500 MeV for the strange 
quark) is used in the quark coalescence model Q, |^ H, Q , 
so the binding energy effect can be neglected. In the 
following, we study the dependence of (f> meson and 
baryon production and their anisotropic flows on the 
strange quark mass. In particular, we use the constituent 
strange quark mass of 500 MeV but same strange quark 
space-time coordinates and momenta at freeze out as be- 
fore. The quark mass effect thus shows up through the 
Lorentz transformation that changes the strange quark 
space-time coordinates and momenta in the rest frame 
of the (j) meson and fl baryon and therefore affects their 
overlap with the quark Wigner phase-space functions in- 
side hadrons. 



Yields and transverse momentum spectra of 1 
and n 



Using a strange quark mass of 500 MeV, we find 
that yields of both (j> mesons and baryons are en- 
hanced significantly compared to previous results with 
a strange quark mass of 199 MeV. For Au-|-Au collisions 
at y/sWN = 200 GeV and 5 = (8) fm, the dN/dy at 
mid-rapidity for (j) mesons with = 0.65 fm changes 



from 8.0 (3.3) to 12.2 (4.8) while that with R^ = 0.47 
fm changes from 4.7 (2.0) to 7.3 (3.0) when the strange 
quark mass changes from 199 MeV to 500 MeV. To keep 
the dN/dy of (j) mesons at mid-rapidity unchanged re- 
quires a reduction of the cj) meson size from i?^ = 0.65 
fm to 0.5 fm or from = 0.47 fm to 0.35 fm when 
the strange quark mass is increased to 500 MeV. For 
baryons, the value of dN/dy at mid-rapidity for 5 = (8) 
fm with i?si = 1.2 fm changes from 0.65 (0.26) to 1.36 
(0.49) when the strange quark mass changes from 199 
MeV to 500 MeV. A somewhat smaller $7 baryon size of 
Rn = 0.9 fm is then needed to keep dN/dy of Q. baryons 
at mid-rapidity unchanged when the strange quark mass 
is increased to 500 MeV. 
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FIG. 16: (Color online) Same as Fig. IHlbut with = 500 
MeV and = 0.5 fm or 0.35 fm and i?n = 0.9 fm. 



We also find that the larger strange quark mass of 
500 MeV leads to a slightly stiffer transverse momen- 
tum spectra for </> mesons and J7 baryons. For R^ = 0.65 
(0.47) fm, the mean transverse momentum (pr) for the 
meson changes from about 0.72 (0.73) and 0.68 (0.70) 
GeV/c for 6 = fm and 6 = 8 fm, respectively, to about 
0.74 (0.76) and 0.69 (0.72) GeV/c for 6 = fm and 6 = 8 
fm, respectively, when the strange quark mass changes 
from 199 MeV to 500 MeV. For baryons, the value 
of (pt) with i?[2 = 1.2 fm changes from about 0.93 and 
0.83 GeV/c for 6 = fm and 6 = 8 fm, respectively, to 
about 0.96 and 0.85 GeV/c for 6 = fm and 6 = 8 fm, 
respectively, when the strange quark mass changes from 
199 MeV to 500 MeV. 

With the new parameters nis — 500 MeV, R,^ = 0.5 
fm or 0.35 fm, and i?si = 0.9 fm fixed from measured 
yields of mesons and baryons, the transverse mo- 
mentum spectra of midrapidity mesons and baryons 
and their ratio in the same collision are shown in Fig. 1161 
As in Fig. |S1 corresponding experimental data are also 
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included in left panels of Fig. UH For = 0.5 (0.35) 
fm and = 500 MeV, the mean (j) meson transverse 
momentum calculated from the 4> meson transverse mo- 
mentum spectrum is about 0.76 (0.77) and 0.72 (0.74) 
GeV/c for 6 = fm and 6 = 8 fm, respectively, which are 
still somewhat smaller than the experimental value of 
about 0.85-1.1 GeV/c For Q baryons with Rn = 

0.9 fm, their mean transverse momentum is about 1.05 
and 0.94 GeV/c for 6 = fm and 5 = 8 fm, respectively. 

The ratio obtained from above new parameters 
is shown in lower panels of Fig. ^| Compared with 
the case of nis = 199 MeV, the ^l/4> ratio increases even 
more from lowpr to intermediate of about 2.5 GeV/c, 
leading to a larger enhancement factor about 14 and 9 
for 6 = fm and 5 = 8 fm, respectively. The r2+/0 at 
PT « 2.5 GeV/c is now about 0.30 and 0.22 for 6 = f m 
and 5 = 8 fm, respectively, which are still significantly 
smaller than that of anti-protons to pions observed in 
experiments but remain comparable to the measured ^l/cjj 
ratio for the larger (j) meson size of Rj, = 0.5 fm. 
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FIG. 17: (Color online) Same as Fig. Clllbut with rris = 500 
MeV and = 0.5 fm or 0.35 fm and Rn = 0.9 fm. 



B. Anisotropic flows of (p and 



strange and antistrange quarks. 



Using the fitted hadron size parameters for = 500 
MeV, i.e., R^ = 0.5 fm or 0.35 fm and Rn = 0.9 fm, 
we show in left panel of Fig. El the pr dependence of 
the anisotropic fiows V2 (solid squares for i?^ = 0.5 fm 
and triangles for i?^ = 0.35 fm) and V4 (open squares 
for i?0 = 0.5 fm and triangles for R^ — 0.35 fm) of 
mid-rapidity 4> mesons produced in Au-|-Au collisions at 
Y^T/VTV = 200 GeV and 6 = 8 fm. Also shown in left panel 
of Fig. El are 1.0w| (solid hne for i?^ = 0.5 fm) and 1.1 
(dashed line for R^ = 0.35 fm). Similar results for 
baryons are shown in right panel of Fig. 1171 Compared 
to those shown in Fig. 1101 with a strange quark mass of 
199 MeV, the hadron elliptic flow is seen to depend only 
weakly on the strange quark mass, with the larger strange 
quark mass of nis = 500 MeV giving a slightly smaller 
V2 at lower px while a slightly larger V2 at higher px- 
Also, the anisotropic flows of 4> mesons and fl baryons 
from the dynamical quark coalescence model with the 
larger strange quark mass of rUs = 500 MeV still satisfy 
the scaling relation va{pt) ~ ^'iCPr) with similar scaling 
coefficients as those for the strange quark mass of 199 
MeV. 

The pT/?^q-dependence of V2 per valence quarks or 
anti-quarks for mid-rapidity 4> mesons and baryons 
produced in Au-|-Au collisions at ^Jsnn — 200 GeV and 
6 = 8 fm, obtained with the large strange quark mass 
of 500 MeV, i?0 = 0.5 fm or 0.35 fm and Rn = 0.9 
fm, is shown in Fig. El together with the elliptic flow 
of mid-rapidity strange and antistrange quarks (s -I- s) 
at freeze-out. It is seen that the V2 of mid-rapidity (p 
mesons and £7 baryons still satisfies approximately the 
valence number scaling in the dynamical quark coales- 
cence model. However, the valence quark number scaled 
V2 remains significantly smaller than the V2 of coalescing 
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Dynamical Quark Coalescence Model 
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n (R =0.9 fm) 




0.0 0.2 0.4 0.6 0.8 
p/n (GeV/c) 



1.0 



FIG. 18: (Color online) Same as Fig. [TTlbut with rUs = 500 
MeV and R4, = 0.5 fm and 0.35 fm as well as Rn = 0.9 fm. 



VII. SUMMARY 

Based on the parton phase-space information obtained 
from a multi-phase transport model within the string 
melting scenario, we have studied the production of (f> 
mesons and fl baryons and their anisotropic fiows in 
Au-f Au collisions at RHIC using a dynamical quark co- 
alescence model, which requires information on the radii 
of (j) meson and f2 baryon. Fixing their radii by fitting 
measured yields of (j) mesons and baryons at midra- 
pidity in central Au-t-Au collisions at -^/stvtv = 200 GeV, 
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we have evaluated their transverse momentum spectra in 
the same colhsion at impact parameter b — and 8 fm 
and also their anisotropic flows in the same collision at 
impact parameter b — 8 fm. Comparing with its value at 
low transverse momenta, we have found that the ratio of 
the yield of Q baryons to that of mesons is enhanced 
significantly at intermediate transverse momenta as ob- 
served in experiments. 

We have further found that the elliptic flows of 
mesons and f2 baryons follow approximately the valence 
quark number scaling. The valence quark number scaled 
elliptic flows of (f> mesons and baryons deviate, how- 
ever, strongly from the underlying V2 of strange and an- 
tistrange quarks. Moreover, we have also studied the 
forth-order anisotropic flow V4 and found that the scaling 
relation of Vi{pT) V2{pt) observed experimentally for 
charged hadrons is satisfied by (j) mesons and 51 baryons 
as well. It will be very interesting to compare these pre- 
dictions with experimental data that are being analyzed. 

In addition, we have studied the dependence of above 
results on the radii of (j) meson and fl baryon as well as 
the strange quark mass. Both the yields and anisotropic 
flows of (f) mesons and f2 baryons are found to be sensitive 
to their radii. For sufficient large radii, the valence quark 
number scaled elliptic flows of (j) mesons and f2 baryons 
approach the V2 of strange and antistrange quarks as in 
the naive momentum-space coalescence model. Also the 
scaling coefficient V4{pt)/v2{pt) is sensitive to hadron 
sizes. Although the strange quark mass was found to af- 
fect significantly the yields of (j) mesons and 57 baryons, it 
does not change much their anisotropic flows. Our results 
thus suggest that in using the quark coalescence model 
to extract the parton dynamics in relativistic heavy-ion 
collisions from hadron observables, it is important to take 
into account the quark structure of hadrons. 

Although results from present study reproduce reason- 
ably the observed transverse momentum dependence of 



the f2/(/) ratio, the transverse momentum spectra of 
mesons and 57 baryons are too soft compared with mea- 
sured ones. This has been attributed to the small current 
quark masses used in the AMPT model. Recent studies 
have shown that the equation of state of the QGP from 
lattice QCD calculations can only be reproduced by par- 
tons with large masses (67l ■ It is thus of great interest 
to improve the AMPT model by using massive partons 
and to study how the resulting parton dynamics affects 
the results obtained in present study. Although we have 
shown in present study that the anisotropic flows of 
mesons and 51 baryons are not sensitive to the strange 
quark mass used in the dynamical coalescence model if 
same strange quark distribution is used, using massive 
partons in the AMPT model not only is expected to af- 
fect the transverse momentum spectra of partons but also 
may influence their anisotropic flow and thus that of pro- 
duced hadrons. Also, it is important to include in future 
studies the contribution from high momentum jets |69l | 
as they are expected to affect the production of mesons 
and 51 baryons as well as their anisotropic flows at inter- 
mediate and high transverse momenta. 
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